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Introduction. The paper presents a theoretical study on binary 
salt ion transport considering the water  dissocia- 
tion/recombination reaction. The work objectives are as fol- 
lows: to build a mathematical model; to develop an algorithm 
for the numerical solution to the boundary value problem cor- 
responding to the mathematical model; to work out the simi- 
larity theory including the transition to a dimensionless form 
using characteristic quantities; to determine a physical mean- 
ing of trivial similarity criteria; to find nontrivial similarity 
criteria; to build and analyze the volt-ampere characteristic 
(VAC). 

Materials and Methods. The theoretical study and numerical 
analysis of the transport of binary salt ions consider the disso- 
ciation/recombination reaction of water. In this case, the heat 
transfer equation and the mathematical model of electrodiffu- 
sion of four types of ions simultaneously (two salt ions, as 
well as H* and OH7ions) in the diffusion layer of elec- 
tromembrane systems with a perfectly selective membrane are 
used. For the first-order differential equations, a singularly 
perturbed boundary-value problem is set. In the equation for 
the electric field, the right side is independent of the intensity. 
In the numerical solution to the digitized system of equations 
by the Newton-Kantorovich method, this causes the stability 
of the method. In this regard, the boundary-value problem is 
reduced for numerical solution: a transition to a system of the 
second-order equations is provided, and the missing boundary 


Beedenue. CraTbaA MOCBALIeHa TeOpeTHYeCKOMy U3yYeHHIO 
TlepeHoca HOHOB OnHApHON COM C y4eTOM peakIMn TMCcouM- 
alu (pekoMOnHalMn) Bob. Lem uccieyqoBaHua: MOcTpoe- 
Hue MaTeMaTHYeCKONH MOeIM; pa3spadoTKa alIropuTMa 4uc- 
JICHHOrO pellleHHA KpaeBOM 3ajjJauH, COOTBeETCTBYIOMIeH MaTe- 
MaTHY4ecKOH MOJeIM; pa3sBuTHe Teopuu Toq00nA, BKIOUAA 
Tlepexoy, K Oe3pa3MepHOMy BUY C MCIOUb30BaHHeM Xapak- 
TePHBIX BeJIMYHH; OlIpesesweHve PU3M4eCKOrO CMBICIAa TPUBHU- 
aJIbHBIX KPHTepveB TOOONA; HaxoxKTeHHe HeETPHBUAIIbHBIX 
KpHTepveB NosoOus; MOCcTpoeHve U aHasIv3 BOJIbTAaMMepHOU 
xapakTepucTuKH (BAX). 

Mamepuaabi u memoooi. IIpu Teopetu4eckom UccleqOBaHun 
WM 4AMCJICHHOM aHasIu3e WepeHoca HOHOB OMHapHOn comm y4u- 
TbIBaeTCA peakINA JMCcouMalMU (pekoMOnHalnN) BorbI. IIpu 
9TOM HCIIOUIb3YIOTCA YpaBHeHHe TeIJIOMNPOBOAHOCTH HM MarTe- 
MaTHuecKad MOJelIb IIEKTPOLMMDy3HU OJHOBPeMeHHO UeTHI- 
peX BHJOB HOHOB (JByX HOHOB COM, a Takxxe HOHOB Ht u 
OH”) B anddy3HOHHOM CyIOe IIEKTPOMeMOpaHHBIX CHCTeM C 
WJleasIbHO CeJIEKTHBHOM MeMOpaHoi;i. 

[na muddepenuMasbHbEIx ypaBHeHHH WepBoro mopsyKa cTa- 
BUTCA CHHTYJIAPHO BO3MYIIeHHad KpaeBasd 3ayaua. B ypaBHe- 
HUM JIA HaNpsKCHHOCTH JICKTPHYeCKOrO MOJIA WpaBad YacTb 
He 3aBHCHT OT HalipsKeHHOocTH. IIpu 4ucneHHOM pellieHuu 
J[UCKpeTH3MpOBaHHOM CHCTeMbI ypaBHeHHi MeTOOM Hbt0To- 
Ha — KaHToposuya 3TO OOycIOBIMBaeT MpoOeMbI yCTOM4n- 
BOCTH MeToya. B cBa3M C STHM KpaeBad 3ajjaua IIpHBONUTCA K 
BUY, YIOOHOMY JIA YMCMeCHHOTO pellleHua: OOecteunBaerca 
llepexoyJ, K CHCTeMe ypaBHeHHM BTOporo MopsyKa, paccuuThI- 
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conditions for the electric field strength are calculated. 
Research Results. A new mathematical model, a numerical 
algorithm to solve a boundary value problem, and software are 
developed. A numerical analysis is carried out, and fundamen- 
tal laws of the transport of salt ions are determined considering 
the dissociation/recombination reaction of water molecules, 
temperature effects, and Joule heating. The VAC is built and 
analyzed. 

Discussion and Conclusions. The transport of binary salt ions 
through a diffusion layer near a cation exchange membrane is 
considered. A mathematical model of this process is proposed. 
It takes into account the temperature effects due to dissocia- 
tion/recombination reactions of water molecules and Joule 
heating in a solution. The basic laws of the transport of salt 
established 


tion/recombination reaction of water molecules and tempera- 


ions are considering the  dissocia- 
ture effects. The temperature effects of the dissocia- 
tion/recombination reaction and the Joule heating in the elec- 
troneutrality region (ENR) are almost imperceptible (with the 
exception of the recombination region, RR). The Joule heating 
in the space-charge region (SCR) is by two orders of magni- 
tude larger than the cooling effect of the water dissociation 
reaction. Upon recombination, approximately the same heat is 
released in the RR as during Joule heating in the expanded 
SCR. However, due to the small size of the RR, the effect of 
this heat is imperceptible. Therefore, we can assume that there 
is only one heat source at the interface in the SCR, which, due 
to its noticeable size, causes a significant increase in tempera- 
ture in the entire diffusion layer. It follows that the emergence 
and development of gravitational convection is possible. Gen- 
eral conclusions, following from the results obtained, open up 
the possibility of intensifying the process of transport of salt 
ions in the electrodialysis machines. 


Keywords: membrane system, ion exchange membrane, space 
charge, extended space-charge 
recombination reaction of water molecules. 


region, dissociation- 
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BalOTCA HeOCTarollMe KpaeBble YCIOBUA JIA HallpsxKeCHHOCTU 
9JIEKTPHYeCKOLO MOA. 

Pe3yibmamol ucciedoeanua. Pa3spadoTaHbl HOBad MaTeMaTH- 
yecKad MOJIeJIb, aJITOPHTM 4HCIICHHOrO pellieHua KpaecBoli 
3alaun, WporpaMMHoe oOecneyeHne. IIpopezeH uncNeHHbIM 
aHasIu3 UW OlpeyieseHbl PyHaMeHTaIbHble 3aKOHOMepHOCTH 
TlepeHoca HOHOB COJIM C y4YeTOM peakIMU AMCcomMalun (pe- 
KOMOMHAaIMH) MOJICKyJI BObI, TEMMepaTypHbIX 39~MeKTOB U 
JIkoysreBa pa3orpesa. IlocrpoeHa u Ipoanasm3upoBaHa BAX. 
O@cyacdenue u 3akiioyeHue. PaccMOTpeH TepeHOC MOHOB 
OuHapHOH comm 4epe3 DHMdy3HOHHbIM CHOU y KaTHOHOOO- 
MeHHOM MeMOpauHpl. IIpeqoxeHa MaTeMaTH4eCKad MOJIeJIb 
Ha3BaHHOrO Wpollecca, KOTOpad Y4NTbIBaeT TeMMepaTypHbIe 
3PbekThI, OOYCNOBIICHHbIC peakKUMAMU JMCCOMMAalHu (peKom- 
OuHalMu) MOseKyJ BOZbI U J>KOyIeBLIM HarpeBOM B pacTBO- 
pe. YcTaHOBJICHbI OCHOBHbI€ 3aKOHOMepHOCTH TepeHoca 
WMOHOB COJIM C yYeTOM peakIIMH MCcolMalun (peKomOnHa- 
I[MH) MOJIeKyJI BOJIbI MH TeMIepaTypHbIX 39PdekToB. Temnepa- 
TYpHble IPekThI OT peakIIMH WHCcoUMalMM (pekKoMOnHalMn) 
u JxKoynepa pa3sorpepa B ODH mpakTnyeckv He3aMeTHBI (HC- 
KIHOUeHHe — OONacTbh pekoMOuHarun, OP). J>koynes HarpeB 
B OOacTH MpoOcTpaHcTBeHHOro 3apsya (OIT3) Ha WBa NopsyKa 
OoNbIe OxakKatomlero osddekta peakiiMu Wuccommalun 
Boybl. IIpu pekomMOuHauuu B OP BblyelaeTcA MpuMepHo 
CTOJIbKO 2Ke Telia, CKOIbKO pu J>KOyIeBOM Harpese B pac- 
mupeHHon OI13. OnHako u3-3a MalIbIX pa3MepoB OP BimmaAHHe 
39TOrO Telia He3aMeTHO. 3HAYHT, MOXKHO CUMTATb, YTO CCTb 
TOJIbKO OJMH MCTOYHUK Telia Ha Mexda3HOu rpannise B OFT3, 
KOTOpbIM OaroyapsA CBOeMY 3AMeTHOMY pa3Mepy OOycIOBIN- 
BaeT 3HAYNTCJIbHOe MOBbILNeCHHe TeEMMepaTypbI BO BCeM JHd)- 
(y3HOHHOM cyIoe. OTcroya cileqyeT, YTO BO3MOX%KHO BO3HHK- 
HOBeHHe MU pa3BHTHe rpaBMTalMOHHOM KOHBeKIIMH. OOmme 
BBIBOJIbI, CJICYIOWMe W3 MOJYYCHHbIX Ppe3yJIbTAaTOB, OTKPbIBa- 
€T BO3MOX%KHOCTb HHTeHCHUKALINUH Wpolecca MepeHoca HOHOB 
COJIM B 9ICKTPOAMasIM3HbIX allnapatTax. 


Ku1r04eBbie CJI0Ba: MeMOpaHHbile CHCTeMbI, HOHOOOMeCHHAaA 
MeMOpaHa, MIPOCTpaHCTBeHHbI 3apay, paciiiu“peHHad OOacTb 
IIPOCTPpaHCTBeHHOTO 3apsyia, WuccouMauua (pekoMOuHalins) 
MOJICKYJI BOJIBI. 


Oopaszeu OA UUumuUpoeanua: MojenupoBaHie UW 4MCJICHHbIM 
aHaiw3 BJIMAHMA peakuMu MccolMaluMu (pe-KOoMOnHaIMN) 
MOJICKYJI BOJ[bI Ha TepeHOC HOHOB com B Auddy3HOHHOM 
cnoe / H. O. Uydsrps [u gp.] // Becruuk Jon. roc. TexH. yH-Ta. 
— 2019. — T.19, Ne3. — C. 268-280. https://doi.org/ 
10.23947/1992-5980-2019-19-3-268-280 


Introduction. Intensive current mode involves using currents several times higher than the limited elec- 


trodiffusion. In this case, secondary (or conjugated) phenomena of concentration polarization occur in the membrane 


systems: 


- spatial electric charge occupies a macroscopic region comparable to the thickness of the diffusion layer; 
- the solution pH changes [1-3], which is explained by the water splitting; 


- microconvective flows are formed, etc. 


Considering impact of the dissociation (recombination) reaction of water molecules is essential to understand- 


ing the transport processes in the electromembrane systems, as some authors believe that the occurrence of new charge 


carriers H* and OH™ can cause a decrease or even disappearance of the space charge, which is the basis for other trans- 


fer mechanisms ( e.g. electroconvection). The dissociation reaction of water molecules is endothermic, and the recom- 


bination reaction 1s exothermic. Centers of the reaction sites are separated in space. This causes an uneven temperature 
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distribution and, therefore, affects all the physical properties of water: dissociation coefficient, heat capacity, etc. In 
addition, an uneven temperature distribution can cause gravitational convection. Thus, the study of temperature effects 
associated with the reactions of dissociation and recombination of water molecules is an urgent task [4]. In this paper, 
we confine ourselves to considering the effect of the uneven temperature distribution on the dissociation coefficient of 
water molecules. 

Materials and Methods 

1. Mathematical model 

1.1. Simultaneous equations 

Consider a one-dimensional stationary transport process for single-charged salt ions in the Nernst diffusion 
layer. We take into account the dissociation (recombination) reactions of water and the temperature effects associated 
with the Joule heating of the solution. Such a case is described through a system of equations [5]: 


=F d dC 
Jt =~ Rr il ei tne (1) 
F d dc 
| RT a7, Dobe T a — D, (2) 
—F a dC 
J3 > RT 343 ax 3 ae (3) 
F dC 
ja= RT pc," . —D, re (4) 
d2 
Te 1 + 3-4), (5) 
0 = AER , i= 1,...,4, (6) 
Rt +E-I+ ee k,7)) =O. (7) 


Here, Rk; = Ro = 0,2, =P, = kaCuo — k,-CaC, = kp(Ky — C34), ky (T) = Kye? -™, i = 1,2 are salt ion 
indices; i = 3 and i = 4 are, respectively, the indices of hydrogen ions H* and hydroxyl OH~; E = —- is the electric 
field strength; C;, j;, Dj are, respectively, concentration, flux, diffusion coefficient of the i-th ion; k,(T) is the constant 
of water dissociation rate; k,. is the constant of the recombination rate of hydrogen ions Ht and hydroxyl OH; k,,is the 
equilibrium constant; ¢€, 1s the dielectric capacitivity; F is the Faraday number; / is the total current density; k is the 
coefficient of thermal conductivity; Po is the density of the solution; c, is specific heat of the solution; G = E - I is the 
density of heat sources associated with the passage of electric current through the solution (Joule heating of the solu- 
tion); B = qk,, Q = B(C; C= <at20) = qk, (C3 = k,,(T)) is the density of heat sources and sinks associated with 
the recombination and dissociation of water molecules. 

The total current density / is equal to: ] = I, + ly. 


Here, J, is partial current density by salt ions; J, is partial current density by water ions. 
From (6)— 8) it follows: j, = const, j2 = const, I, = F - G, — jz) = const, I, = const, but j3 # const, 





j, # const, a? = is thermal diffusivity of the solution. 


Polp 
1.2. Boundary conditions 
Assume that x = 0 corresponds to the depth of the solution, where the electroneutrality condition C, — Cz + 
Cz — C, = 0 is fulfilled, and x = 5 corresponds to the conventional “solution — membrane” interface. Thus, 5™ is the 
thickness of the Nernst diffusion layer. 
In this paper, for definiteness, a cation-selective membrane (CSM) is considered, which for simplicity will be 


assumed to be ideally selective Gy = 0). 





Fig. 1. Diagram of diffusion layer near cation-selective membrane 


1) Boundary conditions in the depth of the solution (x = 0) 
Assume,C,(0) = Cy9, C2(0) = C9, C3(0) = C9, Cy (0) = ae, while Cy9 + C39 = Coo, T(6) = Tp. 
30 
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2) Boundary data at the “solution — membrane” interface (x = 6) 

Suppose that C,(6) = Cy, C3(6) = C3. At this, the sum C,,,, + C3,, is comparable to the tabular value of the 
exchange capacity of the cation-selective membrane. The flux density j,(6) = j,,, of ions OH is determined by the 
catalytic reaction of dissociation of water molecules in a thin surface layer of a cation exchange membrane. 

For temperature, the condition of flow balance is used at the point x = 6 [1]: 


kT = G(8)5 + Q(8)5 + gE 


2. Transition to dimensionless eeu 
2.1. Characteristic quantities and formulas for the transition to a dimensionless form 





As a characteristic concentration of Cg we take the concentration of salt cations C,,. Actual values of Cy, are 
within the range of 10~°M — 1M. As a characteristic temperature, we take the room temperature Ty = 298K. As a char- 
2D1D2 
D4 +D2 





acteristic length, we consider the thickness of the diffusion layer 5, which is of the order of 10~*m. Denote Dy = 


as the electrolyte (NaC/) diffusion coefficient @y = ~~ — as the thermal potential [5]. 


We pass to dimensionless quantities: 
i) =p) = Gp) 88 7 230 P= =e pO = ae —~w pM & ;wW_ Sid _ 











0’ dx’ Inp Inp | Cio ~! DjCio 
T _ a 
7, 7 = 123A, = Iw 
Then, k (T)= we = Cu) gb (7-1) where y) = =bT 5. 
10 


Ae, gt) fe 
10F 252’ F2 Do 
2.2. Transition to dimensionless form in equations 
Focusing on [6], we pass to the dimensionless form in the system of equations (1) - (6): 


Assume that ¢&™ = 


(u) 
ac — plu) -(u) 
oe ee (9) 
(u) 
dc : 
2 = NEM — jt, (10) 
(u) 
dC3 _ p(u) -(u) 
ee (11) 
(u) 
aCy (u) -(u) 
a as (13) 
(u) 
ee eae Oe eee (15) 


We write the total current density in the form of J Se + pe = 1. 


Here, 


j@ = as be _ 1a a pf) = aay ag _ Dj. 


moreover, for a perfectly selective membrane/, jw) — and, consequently, 
[MH = Dp; (Cu) 
S 





Jy 
Under the transition to a dimensionless form, the heat equation (7) takes the form: 
wa (u) . 7) 4 ew (pW _\m),bM(rM-1)) _ 
re ee (CMC? - ye )=0. (16) 
Here, kM = and B™ = aK 8° C10 
DoCioR DoRTo 


2.3. Physical meaning and evaluation of the trivial similarity criteria 

The parameters ec) pO) yO, a Kw, i ) BM, b™ are trivial similarity criteria. Let us determine their 
physical meaning and values at the characteristic input parameters of the problem. 

1. b™ = 15, 

2. The parameter y™ is the square of the ratio of the concentration H* or OH™ ions in a neutral solution 
(pH = 7) to the electrolyte concentration. Its values vary from 107~*° to 1077°, and it can be considered a small param- 
eter. 





2 
3, ew = Aer 9 (=) , where ly = oe is Debye length. Thus, the parameter e™ is the doubled 
10 


C10 F252 %) 
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square of the Debye length ratio I, to the diffusion layer width 5. The values of e™ vary from 10711 to 1077, and this 
parameter can be considered small. 


k k lay | , ; ; 
HG) a Peers, Mlaen ate. (=) . Here, lt, = — . This value is proportional to the width of the re- 
r“10 


F2 Dy  2Ci9F2 Do i 








combination region (RR). Thus, the parameter a™ is the square of the ratio of the Debye length to the width OP I , 
The parameter a“ = 14 and in the framework of this problem, it is a universal constant since it does not depend on the 


input parameters of the problem. So: 7 = Ja) = 3.74, i.e., the width of the OP is approximately 37 % of the width of 


r 
the quasi-equilibrium region of the space charge (SCR, Debye length). 
2.4. Nontrivial similarity criteria, physical meaning and evaluation of quantities 
In the equation (16), the dimensionless coefficients k™ and B™ are trivial similarity criteria. Any relation be- 
tween trivial criteria is called a nontrivial similarity criterion. Below, we express k™ and B™ through other known 
similarity criteria. 
k a*pocp _ a? Poly _ Potp 


2 

. Qa” . . . 

1) Consider k™ = = = . Here, — is the Lewis - Semenov number. We in- 
DoCioR DoCioR Dp CigR CioR Do 














troduce the dimensionless number p = a Then, k™ =- - The number wp indicates the relative specific volumet- 
0Cp 


ric heat capacity, varies from indicates the relative specific volumetric heat capacity, it varies from 1078 to 10, 
and it can be considered as a small parameter. 





qk, 57C q6* 2k-C q 5 5? : 
2) Consider the parameter B™ = Ht? —109 = PA ’rito - 4 Arn -— Here, Arn = + = 22.8 is the 
DoRTo  2RTy Do 2RTo —Po 212 RTo 
r+10 


dimensionless quantity which is a thermal analogue of the Arrhentus criterion. It characterizes the sensitivity of the heat 

; ; 5 lave (0 \" lq\? lq\? 

amount of a chemical reaction to temperature changes. As ~(>) = ~(-4) (=) = = and e=2 (=) a2 (=) : 
2 Ly 2 Ly ld € 5 ly 


then B™ = Arn - = Thus, the heat equation (16) has the following dimensionless form: 




















ar) 
eLe 
3) The calculations show that 
_ kw RR _ kweF*d* R _ kwF?5? 
i C10 Pocp RT9Er Pocp PoCpToEr 
' (uw) la\* kwF2 82 3 os _— ae 
Since es"? = 2 (=) , then N, = ue = 2,5-10°°. Therefore, we obtain a nontrivial similarity criterion, 
oCploEr 
which shows the relation 
yu = née. 
Likewise, 
EU — Nai: 
Here, 
_ R2To Er _ —14 
Ngi = pocpr26? ~ 3,5-10°%, (17) 
7, >= Na 2; (18) 
F*5*k 
Na2 = Renee 
Note that nZ = ng1Nq>.- 
2.5. Dimensionless form of a system of equations 
Thus, a dimensionless system of equations describing the model has the form (index u is omitted for simplici- 
ty): 
dc . dC . dc 
= (iE — jy, = —C2E — jo, oS = C3E — js, 
d , b(T- dC ; 
ED; — jg = a(ye?@-® — C34), = —CyE — jg 


d ., mn dE 
ED, — ja = a(ye’™-Y — €,C,), e— = C,—C, + C3 — Cy. 


2 
eLe — + epEI + padrn(C,C, — ye? TY) = 0. 
The system of equations contains three small parameters €, u and y, connected by relations (17) and (18). The 
parameter € enters the equations singularly, 1.e., as a coefficient in the derivatives, and the parameter y is regular, so the 
system of equations is simultaneously singularly and regularly perturbed. 
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a a a a a a a a 
2.6. Transition to dimensionless form under boundary conditions 
The dimensionless boundary conditions have the form: 


Cy |x=0 = 1, C2\x=0 = 1+ C50, C3|x=0 = C30, C4lx=0 = C49 = aS T(O) = 1, 
Cileaa = (im C3\y=4 = Ca Jalei = Jam: 
, a d 
—k—T(1) = uD,b, —],(1). 
7 ED) WD4bq = Ja) 


Remark 1. Owing to continuity, the following is valid: 


lw = D3j3m — DaJam = Y3J30 — Vajao- 

3. Numerical solution algorithm 

3.1. Transition to a boundary value problem for a system of second-order equations 

Due to the small parameter ¢ in the equations (12), (14) and (15), the boundary-value problem (9) - (16) is sin- 
gularly perturbed; therefore, the system is inconvenient for a numerical solution [7]. In addition, differential equations 
are of the first order, and a boundary-value problem is set for them. In the equation (15), for the electric field strength, 
the right-hand side is independent from FE. This causes stability problems in the numerical solution to the discrete sys- 
tem of equations by the Newton — Kantorovich method [8]. Therefore, the boundary value problem (9) - (16) should be 
reduced to a form convenient for a numerical solution. To do this, we pass to the system of second-order equations: 








ac, _ od d?C, _ -d d?cz3 od _ a(ye?\F-V)_¢3¢4) 
dx2 dx (CE), dx? dx (CLE), dx2 dx (C3E) eD3 ° 
d?Cy _ -d _ afyePT P= C304) PE _ 5 ej, pH 
dx? dx (CyE) cD, : © ane 7 E(C, I C2) hti2t dx dx’ 
a°T sap 6G=1)) = 
eLe 7x2 + euEl + uaArn(C3C, ye ) 0. (19) 


In connection with the transition to systems of the second-order equation, it is necessary to supplement the ini- 
tial boundary conditions with the missing boundary conditions for the function F. From the boundary conditions, it fol- 
lows that at the point x = 0, the electroneutrality condition C, — C, + C3 — C, = 0 1s satisfied. Due to the continuity, 
with some accuracy, in the vicinity of the point x = 0 x = 0, C, —C, + C3 — C, = 0 1s fulfilled, therefore Ey can be 
calculated from the electroneutrality condition. Adding the equations for C;, i = 1, ...,4 and using the electroneutrality 


S— += (¢,(0) - €,(0)). 


a = (Cy — C, + C3 — Cy) |x=1- 


x= 





condition Cig — Cop + C39 — Cyq = 0, we can obtain the relation E|,,-) = 


, d 
From the equation for EF at x = 1, we have: € = 


Research Results 

3.2. Numerical Results 

To analyze the results of the numerical solution, along with the graphs of the desired functions C;, F, j3, j,,T, we 
consider the graphs of the functions p(x) = C, — C, + C3 — Cy, p(x) = C3C, — ye? -. The function p(x) character- 
izes the distribution of charge density, and p(x) is the deviation from the equilibrium of the dissociation (recombina- 


tion) reaction. 
Line Graph:C3*C4-gramma*exp(b*(T-1)) 
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Line Graph: Dependent variable T 
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Fig. 2. Graph of function p(x) = C, — C, + C3 — C, ((0; 0.93]) (a); graph of function p(x) = C3C, — ye? T-)(b); graph of function 
p(x) = C3C, — ye?- in SCR (c); temperature graph (d) 
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From Fig. 2(a), it follows that the diffusion layer is divided into three regions. The first is the region of electro- 
neutrality (ENR), where p(x) = 0. The second is the space charge region (p(x) > 0). The third is an intermediate layer 
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between them where the charge distribution density has a pronounced local maximum. The SCR itself is divided into 
two parts: extended and quasi-equilibrium (it is very small and adjacent to the interface). In the quasi-equilibrium part 
of the SCR, the charge density distribution obeys the exponential law (not shown in Fig. 2). 

A general view of the graph of the function p(x) (Fig. 2, b) shows that the deviation from equilibrium occurs: 

- in the recombination region which is located inside the ENR; 

- in the space charge region (Fig. 2, c). 

In the first case, p(x) > 0, .e., recombination prevails over the dissociation C3C, > ye?-. In the second 
case, on the contrary, p(x) < 0, and dissociation prevails over recombination in the SCR, as well as in the intermediate 
layer. The temperature (Fig. 2, d) increases linearly by 0.0028 in a dimensionless form (or 0,0028 - 298 = 0.83 de- 
grees). The increase depends on the total current density or potential jump in the diffusion layer and can be up to tens of 


degrees. Obviously, all terms of the equation (19) are small compared to the first one, that 1s, the second derivative is 
ae ., o ue ; a 
= =~ 0. From this, it follows that the temperature is linear and is determined by the boundary conditions. 
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Fig. 3. General view of graphs of functions C,, Cz, C3, C4 (a); function graphs C, Cz near zero (b); function graphs C3, C4 
near zero (c); function graph C,, Cz, C3, C4 in space charge region (d) 


As follows from Fig. 3 a, the sodium concentration decreases linearly in the ENR and becomes small in the 
SCR, but it is higher than the concentration of the remaining ions (Fig. 3 d). The center of the RR is the intersection 
point of the graphs of concentration of Na* u Cl and, respectively, of OH~ and H*. To the right of the center, the 
concentration of Cl~ decreases faster than the concentration of Na*. The deficit of negatively charged ions resulting 
from this is compensated by the rapid growth of hydroxyl ions, and the condition of electroneutrality is preserved (Fig. 


3a). 
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d) 
Fig. 4. General view of graph of FE function with exception of phase boundary (a); graph of F function near RR (4); graph of FE func- 


tion near phase boundary (c), volt-ampere characteristic (d) 
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From the data in Fig. 4, it follows that in ENR, the electric field strength increases rapidly tens of times (Fig. 
4, a). A particularly high rate of the electric field strength is observed in the SCR upon reaching values of the order of 
1 
Je 
carriers change: before RR, the charge is determined by Na*, Cl~, H* ions, and after — by Na*, Cl-, OH “ions. The 
appearance of a “step” owes to the fact that the diffusion coefficient of hydrogen is approximately two times greater 
than the diffusion coefficient of hydroxyl. The VAC increases without limit with a growth of the potential jump in the 
over-extreme mode J, > 2 (Fig. 4, d). This growth tends to be much slower than the experimentally observed one, 
which means that along with the dissociation (recombination) reaction of water, there is another mechanism of super- 
limit transfer. It 1s generally accepted that electroconvection is such a mechanism [9-11]. In this regard, a problem of 
evaluation and interaction of the dissociation (recombination) reaction and electroconvection occurs. 


In the RR there is an internal boundary layer of F electric field strength (Fig. 4, b). This is due to the fact that charge 
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b) 
Fig. 5. Flux plot of hydrogen j, (a) and hydroxylj, (b) ions 


Fig. 5 shows that the flow of hydrogen ions is positive. It is constant to the left of the RR and is practically 
equal to zero to the right of the RR till the SCR. In the SCR, the flow of hydrogen ions increases linearly due to the in- 
tense dissociation of water. The flow of hydroxyl ions to the left of the RR 1s practically zero, and to the right, it is 
negative and constant until the SCR. In the SCR, the flow of hydroxyl ions, as well as the flow of hydrogen ions, in- 
creases linearly. In a quasi-equilibrium SCR, the hydroxyl ion flux changes sharply and satisfies the boundary condition 
jaQ) = ja(m). The recombination region is an internal boundary layer for the flows of hydrogen and hydroxyl — they 
change sharply here. 

Discussion and Conclusions. In this paper, the transport of binary salt ions through a diffusion layer near a 
cation exchange membrane is considered. A mathematical model of the process is proposed, which takes into account 
the temperature effects caused by the dissociation (recombination) of water molecules and the Joule heating in solution. 
This model consists of 9 ordinary differential equations with corresponding boundary conditions. An algorithm for the 
numerical solution to a boundary value problem is developed, and its computational investigation is presented. The 
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basic laws of the transfer of salt ions are defined considering the dissociation (recombination) reaction of water mole- 
cules and temperature effects. 

It is shown that the VAC grows without limit with a growth of the potential jump in the over-extreme mode, 
and is much slower than in experiments. Therefore, along with the dissociation (recombination) reaction of water, there 
are other mechanisms of super-limit transfer (for example, electroconvection and gravitational convection). Thus, a 
challenge arises of assessing the interaction of dissociation (recombination) reactions, electroconvection, and gravita- 
tional convection. 

Temperature effects from the dissociation (recombination) reaction and the Joule heating in the ENR are al- 
most imperceptible (with the exception of RR). The Joule heating in the SCR is two orders of magnitude higher than the 
cooling effect of the water dissociation reaction. Upon the recombination, approximately the same heat is released in the 
RR as during the Joule heating in the expanded SCR. However, due to a small size of the RR, the effect of this heat is 
imperceptible. Thus, we can assume that there is only one source of heat at the interface in the SCR, which, due to its 
noticeable size, causes a significant increase in temperature throughout the diffusion layer. It follows that the emergence 
and development of gravitational convection 1s possible. 

In the paper, a solution of sodium chloride is considered; however, the study is valid for a solution of any bina- 
ry salt. 

General conclusions following from the results obtained enable to intensify the process of salt ion transport in 
the electrodialysis machines. 
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